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ABSTRACT 


The  development  of  a  high  average  power  FEL  for  military  applications,  whether 
shipboard  or  not,  represents  a  significant  advancement  in  technology  over  present 
weapon  systems  design.  The  FEL  has  significant  advantages  over  conventional 
kinetic  systems  and  other  classical  high-energy  laser  systems.  The  rapid  response, 
wavelength  tunability,  and  infinite  magazine  make  the  FEL  a  highly  desirable 
shipboard  weapon  system. 

The  initial  part  of  this  thesis  examines  the  advantages  of  a  FEL  over  a 
conventional  kinetic  weapon.  Section  II  explores  the  atmospheric  phenomena  that 
affect  the  propagation  of  a  laser  beam  enroute  to  its  target.  Section  III  presents  the 
Boeing  FEL  proposal  followed  by  the  theory  of  the  FEL.  Lastly,  in  Sections  V  and  VI, 
simulations  are  conducted  to  analyze  the  FEL’s  feasibility. 
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I.  INTRODUCTION 


During  the  Cold  War,  the  "Super  Powers"  (U.S.  &  U.S.S.R.)  amassed 
enormous  stockpiles  of  strategic  and  tactical  weapons.  Versions  of  these 
weapons  were  distributed  to  third  world  allies  by  both  Super  Powers.  As  the 
cold  war  dragged  on,  developed  countries  such  as  China  and  France 
dramaticaiiy  increased  their  missile  design  capability.  Accordingly,  these 
nations  began  to  sell  highly  advanced  missile  systems  to  the  highest  bidders. 
With  the  demise  of  the  cold  war,  the  single  Soviet  adversary  no  longer  existed. 
What  resulted  was  the  creation  of  many  potential  enemies  throughout  the  world. 
At  present,  any  nation  or  organization  with  resolve  and  financial  means  can 
obtain  and  use  modern  offensive  missile  systems.  Therefore,  we  as  a  nation, 
as  well  as  our  democratic  allies,  are  susceptible  to  lurking  threats  from  a 
number  of  conceivable  directions.  Political,  religious  or  terrorist  groups  of 
various  sizes,  as  well  as  developing  countries  struggling  to  take  a  stand  are  all 
capable  of  attack. 

During  Operations  Desert  Shield  and  Desert  Storm,  a  major  concern  of  the 
coalition  forces  was  the  possible  employment  of  Chinese-made  Silk  Worm  Anti- 
Ship  Missiles  (ASM)  by  Iran.  Fortunately,  no  ASMs  were  fired  at  the  coalition 
forces.  In  1987  however,  during  the  Iran-lraq  War,  the  Iraqis  inflicted  significant 
structural  damage  on  the  U.S.S.  STARK  [1]  when  they  successfully  fired  a 
French  made  Exocet  ASM.  In  a  similar  missile  incident  during  the  Falklands 
War,  the  Argentinian  Navy  inflicted  catastrophic  damage  on  a  British  Frigate 
from  an  air  launched  Exocet  missile  [1].  The  world  learned  two  lessons  from 
both  these  incidents:  first,  the  potential  for  destruction  from  an  ASM  was 
greater  than  anticipated;  second,  vessels  at  sea  are  vulnerable  to  darnage 
because  present  Anti-Ship  Missile  Defense  (ASMD)  technology  is  far  inferior  to 
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ASM  technology.  Moreover,  we  came  to  the  alarming  realization  that  the 
defensive  systems  utilized  by  the  British  and  U.S  Forces  failed  to  protect  the 
warships  from  damage.  The  United  States  relied  heavily  on  the  Close  In 
Weapons  System  (CIWS)  -  essentially  a  20mm  Gatling  Gun  -  and  the  British 
relied  on  the  Penguin  Gatling  Gun. 

ASM  technology  continues  to  dramatically  improve  while  the  predominately 
1960’s  technology  of  present  Anti-Ship  Missile  Defense  System  (ASMDS)  have 
only  improved  in  the  detection  aspects  of  the  system.  As  an  example,  the 
French  will  soon  deploy  the  Aerospatiale/MBB  Supersonic  Anti-Ship  Missile 
(ANS)  [2]  to  replace  its  Exocet  missile  stockpile.  The  ANS  will  be  a  long  range 
(  =180  km)  sea  skimming  missile  capable  of  10  g  terminal  deceptive 
maneuvers  (  1  g  is  9.8m /sec  ^  )  at  speeds  in  excess  of  mach  2.  As  a  rule  of 
thumb,  ASMD  systems  need  to  maneuver  at  least  3  times  more  rapidly  than  the 
offensive  ASM  [3].  For  a  mach  2  missile  with  10  g  maneuverability  this  would 
require  an  ASMD  System  that  could  maneuver  at  30  g.  If  the  technology 
existed  to  make  an  ASMD  system  maneuver  at  30  g’s,  in  time  it  would  be 
incorporated  into  the  ASM  design.  The  endless  advancement  cycle  would 
always  benefit  the  offensive  weapons  system.  Although  the  two 
aforementioned  examples  focus  on  the  failure  of  mechanical  defensive  weapons 
technology  to  guard  against  a  potential  supersonic  ASM,  the  technology  exists 
to  develop  a  high  energy  laser  system  to  serve  as  a  "speed-of-light"  Anti-Ship 
Missile  Defense  System.  This  laser  system  would  terminate  the  endless 
advancement  cycle. 

The  speed-of-light  maneuverability  of  a  Shipboard  Free-Electron  Laser 
(SFEL)  is  just  one  advantage  over  a  kinetic  energy  defensive  weapon  system 
such  as  a  missile.  Other  characteristics  that  make  a  SFEL  a  superior  ASMDS 
are  its  relatively  endless  magazine  supply,  its  tunable  wavelength  and  its  high 
power  potential.  Section  11  of  this  thesis  will  describe  the  amount  of  localized 
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enGrgy  TGCjuirGd  to  dGStroy  a  missilG  and  will  follow  thG  propagation  path  of  thG 
lasGf  boam  back  to  thG  SFEL  to  dGtGrminG  thG  bGam  rGquirGnriGnts  at  thG  lasGr. 
NumGhcal  simulation  is  usGd  to  dsscribs  ths  offsets  of  diffraction,  turbulence, 
and  thermal  blooming  on  laser  beam  propagation.  Section  III  will  describe 
Boeing  Corporation’s  1  Kilowatt  Average  Power  PEL  and  its  potential  for  use  as 
a  SFEL.  In  Section  IV  the  basic  principles  of  the  FEL  are  described.  Lastly, 
simulations  and  basic  FEL  theory  are  used  to  evaluate  the  performance  of  the 
Boeing  FEL  in  Section  V. 
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II.  LASER  BEAM  PROPAGATION  THROUGH 
THE  ATMOSPHERE 


A.  DESCRIPTION  OF  THE  ATMOSPHERE 

The  transmission  of  laser  radiation  through  the  atmosphere  results  in 
several  processes  that  alter  the  nature  of  the  beam.  At  low  average  power 
where  the  interactions  are  linear,  the  effects  of  scattering  and  absorption  are 
dominant.  At  high  average  power,  additional  non-linear  processes  can  occur. 
The  performance  of  a  high  average  power  laser  operating  at  sea-level  can  be 
severely  degraded  by  these  non-linear  effects. 

As  a  model  for  development,  our  firing  platform,  the  SFEL,  and  Exocet 
target  are  both  at  sea-level  and  approximately  14  kilometers  apart.  The  target 
is  in  the  Line-of-Sight  of  the  SFEL  during  the  entire  optical  wave  propagation. 
In  the  absence  of  an  atmosphere  the  power  of  the  laser  beam  would  be 
constant  in  time.  The  atmosphere,  however,  produces  spatial  and  temporal 
fluctuations  along  the  propagation  path,  resulting  in  changes  in  the  intensity  and 
power  of  the  beam.  This  is  due  to  random  variations  in  the  atmospheric 
refractive  index  along  the  propagation  path.  The  variations  in  the  atmospheric 
refractive  index  primarily  result  from  small  temperature  variations  along  the 
propagation  path.  To  better  understand  the  effect  of  the  refractive  index 
fluctuations,  consider  the  region  of  high  or  low  refractive  index  as  a  lens.  In  the 
following  section,  I  will  develop  the  propagating  wave  using  Maxwell’s  equations 
and  describe  the  refractive  index  variations  as  it  applies  to  the  wave  equation. 

The  turbulent  atmosphere  can  then  be  modeled  as  a  large  number  of 
random  lenses,  varying  in  size  and  shape  along  the  propagation  path  [4].  The 
turbulence  will  cause  no  change  in  the  magnitude  of  the  electromagnetic  wave, 
but  a  phase  change  results  from  the  change  in  velocity  of  the  area  of  the  wave 
that  traverses  the  lens.  The  resulting  phase  variations  are  not  uniform 
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throughout  all  points  perpendicular  to  the  direction  of  propagation.  These 
varying  perturbations  may  cause  focusing  or  defocusing  effects  and  local 
deviations  in  the  direction  of  wave  propagation.  Eventually,  through 
interference  and  diffraction,  there  are  intensity  fluctuations  and  loss  of  intensity 
at  the  missile  surface.  While  the  individual  effect  of  a  small,  single  refractive 
index  change  over  a  small  distance  is  minute,  the  cumulative  effect  of  many 
refractive  index  variations  can  be  significant. 


B.  DIFFRACTION 


To  analyze  the  propagating  wave,  it  is  essential  to  begin  with  Maxwell’s 
Wave  equation  [5].  An  electromagnetic  wave  will  propagate  in  the  absence  of  a 
current  density  according  to 
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{t,t)  =  0 


(1) 


The  first  term  is  the  Laplacian  +  dy  +  3|,  t  is  the  time,  c  is  the  speed  of 

light,  /  is  the  vector  position,  and  A  is  the  optical  vector  potential.  The  vector 
potential  is  complex,  containing  both  phase  and  amplitude  information.  The 
optical  electric  and  magnetic  fields  can  be  derived  from  the  vector  potential  by 
the  following  relationships: 

S  =  -  “  ,  (2) 

C  at 


and 

S  =  ^xA  .  (3) 

The  optical  beam  propagates  in  the  z  direction  along  the  range,  R,  to  the 
target.  The  beam  is  assumed  to  be  composed  of  a  single  frequency  and 
circularly  polarized.  The  vector  potential  describing  such  a  beam  is 
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where  k  =  27cc/to  is  the  wave  number  of  the  beam,  and  ^  is  the  phase 
describing  a  plane  wave,  ^  =  kz  —  (ot  +  ^.  In  this  expression,  co  is  the 
angular  frequency  and  (|)  is  the  optical  phase.  Since  the  beam  is  composed  of 
a  single  frequency,  the  variation  in  the  electric  field  amplitude  and  phase  will  be 
small  over  the  distance  of  a  wavelength,  X  =  2n/k.  Thus,  BE  /  dz  kE,  and 
3(|)  /  az  «:  k^.  The  variation  in  the  field  amplitude  and  phase  in  time  over  a 
period  will  also  be  very  small,  so  that  BE  I  Bt  <  coE ,  and  B(^  /  Bt  ^  (ocf). 
Thus,  the  wave  equation  contains  second-order  terms  that  will  be  small 
compared  to  first  or  zeroth-order  terms.  In  order  to  simplify  the  analysis,  these 
second-order  terms  are  dropped.  Substituting  in  the  vector  potential,  the  wave 
equation  then  becomes  [3], 


+  2ik^  +  — ^exp(/(t))  =  0 
ay2  c  Btjj 


The  last  partial  derivatives  in  z  and  t,  describe  the  propagation  of  a  plane 
wave.  The  first  partial  derivatives  in  x  and  y  describe  how  the  beam  changes 
due  to  diffraction  as  it  propagates. 


Another  simplification  is  to  make  the  terms  dimensionless  for  numerical 
analysis.  A  dimensionless  time  is  constructed  from  %  =  ct/R  where  x  varies 
from  0  to  1  over  the  entire  range  of  travel.  It  is  convenient  to  follow  the 
wavefront  with  Z  =  z  -  ct.  With  this  coordinate  change. 


_3_  _  ^_B_  ^_B_  ^  A. 
Bz  ~  Bz  BZ  ^  Bz  9x  BZ 


(6) 


using  B^Z  =  1  (  where  B^  =  B  I  Bz  )  and  B^i  =  0.  Also, 

lA  =  +  lilA  =  ±1 _ A  (7) 

c  Bt  ~  c  Bt  BZ  c  Bt  Bx  R  Bx  BZ 
using  BfZ  =  -  c  (  where  Bt  =  B/Bt  )  and  BfX  =  c  /  R.  The  plane-wave 
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operator  in  the  wave  equation  (5),  now  becomes 


9  1  _  9  9  _  1  9 

'cYt  "  YY  ~  dZ  ~  R  dx 

The  wave  equation  becomes 


dx^  dy^ 


„ 1  9 
+  2//c— — 
R  9t: 


^exp(/(t))  =  0 


(8) 


(9) 


and  is  known  as  the  parabolic  wave  equation.  The  parabolic  wave  equation 
relates  the  transverse  size  of  the  beam  to  its  propagation  distance  R,  and  can 
be  rewritten  as 


^Vf  + 
2k  -‘■ 


9t 


^exp(/(l))  =  0 


(10) 


where  the  Vf  operator  is  the  transverse  Laplacian,  Vf  =  9^  +  9w.  If  the 


transverse  variables,  x  and  y,  are  made  dimensionless  through  x  x^k/2R , 
and  y  y^k/2R ,  Equation  (10)  becomes 


f 


(11) 


where  the  complex,  dimensionless,  optical  field  is  a  =  la  I  exp(/(t)).  Since  (11) 
is  linear  in  the  electric  field,  E  oc  a,  a  simple  scale  factor  can  make  the  field 
dimensionless  and  normalized  to  unity.  Equation  (11)  describes  how  the 
dimensionless  complex  field,  a{x,y,x),  evolves  in  the  dimensionless  time,  t,  and 
the  dimensionless  transverse  coordinates  x  and  y. 

The  complex  field  may  also  be  written  as  the  sum  of  its  real  and 
imaginary  components  a  =  a^.  +  /a,  =  la  I  exp{i(^).  Equation  (11)  can  thus 
be  separated  into  real  and  imaginary  parts. 


5a,  1  2 

(12) 

9t  “  4^^^^  • 

(13) 
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Equations  (12)  and  (13)  describe  the  motion  of  the  laser  beam  as  it 
propagates  to  the  target  and  include  diffraction  of  the  optical  wave.  The 
complex  field  a  evolves  in  time  x  along  the  range  R.  As  x  =  0  — >  1 ,  the 
wavefront  evolves  along  the  range  to  the  target.  A  wavefront  of  large  size, 
Ax,  Ay  ^  a  plane  wave,  will  not  be  affected  by  diffraction  since  the 
divergence  of  a  becomes  small.  However,  if  the  structure  size  is  smaller,  it  will 
be  affected  by  diffraction.  Specifically,  if  Ax  =  Ay  characterizes  the  spatial 
extent  of  the  beam  in  cross  section  and  Ax  =  1  over  the  full  range  of 
propagation,  then  diffraction  is  important  when  Ax  /  (4Ax^)  =  1.  In  summary,  a 
large  wave  front  with  structure  size  Ax,  Ay  »  1  is  not  seriously  affected  by 
diffraction  compared  to  a  small  wavefront  with  structure  size  Ax,  Ay<  1. 
Equations  (12)  and  (13)  are  linear,  complex  equations  with  symmetry  in  the  x 
and  y  directions. 

Equations  (12)  and  (13)  can  be  solved  numerically  and  the  complex  field 
determined  for  any  x  and  y  as  the  beam  traverses  x.  The  propagation 
simulation  code  [3]  uses  two-dimensional  (2-D)  arrays  for  the  real  (ar[x][y])  and 
imaginary  (a/[x][y])  parts  of  the  field.  The  program  then  steps  these  arrays 
along  the  dimensionless  time,  x  =  0  to  1 .  At  each  step,  the  program  calculates 
the  new  part  of  the  field,  real  and  imaginary,  using  the  previous  parts  of  the 
field.  The  Laplacian  is  calculated  by  determining  the  field  difference  between 
the  X  and  y  points  prior  to  and  after  each  time  step.  Every  2-D  array  has 
square  dimensions,  and  each  address  assignment  contains  a  floating  point 
number.  The  mathematical  derivation  of  this  process  is  as  follows: 

a[x][y]  =  ar[x][y]  +  /a/[x][y]  .  (14) 

In  this  derivation  process  a  refers  to  either  ar  or  ai.  Evaluating  to  the  right 
gives 


da 

dx 


a[x  +  1][y]  -  a[x][y] 
dx 


(15) 


and  evaluating  to  the  left  gives 


^  ^  a\x]\v]  -  a[x-1][y] 
dx  dx 


(16) 


To  evaluate  d^a/dx^,  it  is  important  that  both  right  and  left  partial  derivatives  be 
used  to  maintain  the  symmetry  of  Equation  (11)  and  to  balance  between 
adjacent  array  locations.  Thus, 

,  (17) 
(18) 
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'da^ 

1 

a[x  +  1][y]  -  a[x][y] 

a[x][y]  -  a[x  -  1][y] 

3x 

dx 

dx 

dx 

dx 

3fa  ^  a\x  +  11M  +  a[x  -  DIH  -  2a[x][y] 
dx^  dx^ 


The  y  component  equations  are  similarly  derived.  Equation  (11)  can  be  written 
numerically  as 


dar[x]\y]  = 


32  ^  32 


a/[x][y] 


(19) 


and 


da/[x][y]  = 


dx  3^ 
4  [3x2 


ar[x][y] 


(20) 


Where  dar  and  da/  are  respectively,  the  dimensionless  time  derivative  of  the 
real  and  imaginary  parts  of  the  field.  The  numerical  coefficient  C  in  the 
algorithm  is  dx/{4dx^).  It  serves  as  the  propagation  front  diffraction  operator, 
where  the  time-step  is  dx  =  1/nf,  and  nt  is  the  total  number  of  time  steps  over 
the  range  0^1.  Additionally,  the  finite  grid  size  is  dx  =  dy  =  W/nx,  which  is 
the  quotient  of  the  transverse  window  width  (W)  and  the  number  of  bin 
elements  (nx)  in  the  x  and  y  directions.  Substituting  Equation  (18)  and  its  y 
counterpart  into  Equations  (19)  and  (20)  yields  a  numerical  solution  to  the  wave 
equation.  Furthermore,  as  explained  in  the  analytical  derivation,  a  wavefront 
will  be  affected  by  diffraction  if  C  =  1 . 
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An  analytical  solution  to  the  optical  wave  equation  gives  the  following 
relation  between  the  dimensionless  Rayleigh  length  Zq  [3],  the  mode  waist 
radius  Wq,  and  the  wavelength  X: 


Zq  = 


(21) 


The  Rayleigh  length  gives  the  characteristic  distance  over  which  the  beam  waist 
will  double  in  transverse  area  along  its  propagation  path  [6]. 

Figure  1  describes  a  typical  Gaussian  optical  mode  in  a  resonator  with 
spherical  mirrors.  This  is  the  output  of  the  simulation  code  for  Zq  =  0.5  and  the 
location  of  the  mode  waist  in  the  middle  of  the  range  =  0.5.  Where  x^  is  the 
dimensionless  time  which  specifies  the  mode  waist  location.  At  the  laser 
beam  spreads  by  diffraction  as  it  moves  beyond  The  top  left  window 
depicts  the  side  view  of  a  beam  profile  with  the  peak  dimensionless  optical  field 
shown  in  brackets.  The  intensity  scale  of  the  windows  is  such  that  regions 
where  no  light  is  present  are  dark  and  the  bright  regions  indicate  where  light 
has  propagated.  The  contour  lines  show  where  the  optical  field  amplitude  is 
50%  of  the  maximum.  The  vertical  axis  is  the  dimensionless  x  direction  and  the 
horizontal  axis  is  the  dimensionless  time  variable  x  =  0  ^  1 .  The  wavefront 
initially  focuses  at  x  =  0.5  and  then  it  spreads.  As  the  beam  passes  the  mode 
waist  location,  it  spreads  along  x.  The  top  right  window  shows  the  peak 
dimensionless  optical  amplitude  at  the  end  of  the  simulation  x  =  1 .  The  bottom 
window  shows  the  peak  dimensionless  optical  field  from  a  top  view  perspective. 
Observe  how  the  optical  field  is  symmetric  about  the  x  and  y  directions  from 
both  the  top  and  side  views.  This  results  from  having  no  inhomogeneities  along 
the  propagation  path.  The  data  on  the  right  lists  the  value  of  the  numerical 
coefficient  C,  the  number  of  time  steps  nt,  the  number  of  bin  elements  nx,  and 
the  corresponding  error  in  the  calculation.  The  error  computation  is  given  by 
the  ratio  of  the  change  in  energy  to  the  energy  dElE  over  consecutive  bins. 
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Specifically,  this  relationship  is  given  by: 


E  =  . 


f  =  ,  (23) 

tz  Cq 

where  Ef  and  Eq  are  respectively,  the  final  and  initial  energies.  Below  the  error 
calculation  is  an  optical  field  intensity  bar  scale.  The  scale  shows  the 
magnitude  of  the  dimensionless  optical  field  amplitude. 


**  Laser  Beam  Propagation 

T  =0.5  z^=0. 5 

w  o 


top(0.99) 


C=0. 01562 
nt=1024,  nx=64 
dE/E=0.0019 


a  (max) 


Figure  1 .  Propagation  path  with  no  inhomogeneities. 


An  alternate  method  of  solving  Equation  (11)  is  by  implementing  a  Fast 
Fourier  Transform  (FFT)  [6].  The  first  step  in  this  process  is  to  show  that 

'  ixvf 

a(t)  =  exp  a(0)  ,  (24) 
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is  a  solution  to  Equation  (11).  Thus,  substituting  Equation  (24)  into  (11)  yields: 
da 


dx 


'  ixvf 

/Vf 

r 

/xVf 

1  4  JJ 

.=a(0)A 

exp 

-L 

4 

_  X 

4 

J 

a  (O)exp 

;vf 


a(t)  ,  (25) 


which  is  Equation  (11).  Therefore,  Equation  (24)  is  a  solution  to  the 
dimensionless,  parabolic  wave  equation.  Translating  Equation  (24)  by  a 
dimensionless  time  step  Ax  results  in: 


a(x  +  Ax)  =  exp 


/AxV^ 

”4 

V  ✓ 


a(x) 


(26) 


An  exact  numerical  solution  can  then  be  derived  by  applying  a  FFT  to  Equation 
(24)  and  multiplying  it  by  the  diffraction  operator,  and  then  taking  the  inverse 
FFT.  Without  tracing  through  the  mathematics,  the  procedure  is: 


a  (x  +  Ax)  =  FFT  ^ 


exp 


-/Ax/c^ 

4 


FFT{a{x)) 


(27) 


Since  this  solution  is  exact,  and  Ax  need  not  be  infinitesimal,  no  numerical 
errors  are  introduced  into  the  propagation  of  the  wavefront. 

Figure  2  illustrates  a  Gaussian  optical  mode  in  a  resonator  with  spherical 
mirrors  utilizing  the  FFT  technique.  This  is  the  output  for  Zq  =  0.2  and  the 
mode  waist  is  located  at  the  end  of  the  range  x^,  =  1 .0.  As  expected  no  errors 
propagated  through  the  calculation.  This  can  be  seen  from  dE/E  =  0.  In  this 
simulation  the  beam  focuses  at  the  end  of  the  propagation  path.  The  diffraction 
operator  forces  the  beam  to  focus  at  the  beam  waist  location  x^^.  This  is  of 
particular  interest  to  a  SFEL  since  the  maximum  field  intensity  should  be 
concentrated  on  the  target.  The  following  sections  explore  how  the  atmosphere 
affects  beam  quality. 
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C=0. 01562 
nt=1024,  iix=64 
dE/E=0.0000 


Figure  2.  Propagation  simulation  incorporating  FFT  technique. 


C.  MODELING  A  LENS 

As  a  laser  beam  travels  through  a  volume  element  in  the  atmosphere,  or 
in  any  medium,  its  propagation  is  affected  by  the  medium’s  index  of  refraction 
n  =  c/v.  Where  c  is  the  speed  of  light  in  free  space  and  v  is  the  speed  of  light 
in  the  medium.  The  non-uniform  density  of  the  volume  element  results  in  a 
changing  index  of  refraction  across  the  beam  cross-section.  The  varying  index 
of  refraction  along  the  transverse  axis  of  the  beam  is  similar  to  the  change  in 
the  refractive  power  across  a  lens  depending  on  the  thickness  of  the  lens  [3].  A 
lens  will  have  varying  affects  on  distinct  segments  of  the  beam.  The  overall 
effect  is  to  induce  a  phase  shift  across  the  wavefront  which  is  dependent  on  the 
thickness  of  the  lens  it  passes  through.  Specifically,  this  phase  shift  can  be 
represented  by  the  following  relation: 

3out  =  a/nexp(/A(())  .  (28) 
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Where  3^  is  the  complex  field  entering  a  path  segment  of  length  Az,  A<t)  is  the 
phase  shift  resulting  from  the  Az,  and  3^^,  is  the  exiting  complex  field. 
Comparing  a  wavelength  of  a  ray  which  passes  through  the  center  of  the  lens 
with  one  which  passes  through  a  portion  of  the  lens  a  distance  r  from  the 
center,  a  phase  shift  will  be  introduced  between  the  two  rays  according  to  the 
relation: 

A(|,  =  .  (29) 

Where  Acj)  is  the  dimensionless  phase  shift,  and  is  the  radius  of  curvature  of 
the  lens.  The  propagation  code  uses  the  operator  in  Equation  (28)  to  model  the 
effect  of  a  lens  on  the  propagation  of  the  beam. 

As  an  example.  Figure  3  illustrates  placing  a  lens  at  t  =  0.3  with  a 
strength  of  =  3.  Where  is  the  phase  delay  (or  advance)  introduced  by 
the  presence  of  a  lens.  The  portion  of  the  beam  passing  through  the  center  of 
the  lens  is  impeded  more  than  at  the  edges.  The  effect  of  the  lens  is  to  cause 
strong  focusing  of  the  beam  at  x  =  0.6.  At  t  =  1,  the  beam  is  beyond  its  focal 
point  and  spreads  out  again.  Additionally,  this  simulation  incorporates  a  Fast 
Fourier  Transform  technique  to  propagate  the  wavefront.  Figure  4  simulates 
having  the  lens  in  the  same  position  but  with  a  lens  strength  of  d^  =  -3.  In  this 
case  the  portion  of  the  beam  passing  through  the  center  of  the  lens  is  impeded 
less  than  the  edges  and  the  lens  defocuses  the  beam. 

D.  MODELING  TURBULENCE  WITH  LENSES 

Consider  the  SFEL  model  as  the  transmitter  and  a  stationary  platform  as 
the  target  within  the  line-of-sight  (LOS)  of  the  SFEL.  When  the  SFEL  is  fired  in 
a  continuous  wave  mode  at  the  platform  target,  the  intensity  of  the  intervening 
turbulence  the  beam  wave  encounters  as  it  propagates  towards  the  target 
varies  in  a  diurnal  cycle  [4].  Atmospheric  turbulence  effects  are  smaller  at  night 
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top(2«86) 


C=0. 01562 
nt=1024,  nx=64 
dE/E=:0.0000 


I  a  (max) 


Figure  3.  Propagation  path  with  a  converging  lens. 


and  in  the  early  morning  hours  than  during  the  mid-day  hours  when  solar 
heating  of  the  earth  becomes  significant  enough  to  initiate  convective  heating 
currents  [7].  The  convective  heating  current  or  instability  generates  temperature 
fluctuations  as  a  function  of  height.  That  is,  the  sun-heated  surface  materials 
warm  a  thin  layer  of  the  air  above  it.  Air  parcels  from  the  layer,  displaced 
upward,  are  warmer,  less  dense,  and  more  buoyant  than  the  ambient  air.  The 
mixing  of  these  hot  air  parcels  with  the  cool  descending  air  parcels  produces 
isobaric  temperature  irregularities  which  cause  fluctuations  in  the  index  of 
refraction.  The  measurement  of  these  temperature  fluctuations  provides  a 
measure  of  local  turbulence  [8].  If  one  considers  turbulence  as  varying  indices 
of  refraction  through  a  small,  isobaric,  volume  element,  then  turbulence  can  be 
modeled  as  a  series  of  random  lenses  with  varying  indices  of  refraction. 


Figure  4.  Propagation  path  with  a  diverging  lens. 


The  simulation  code  uses  a  series  of  random  lens’  superimposed  on  each 
other  to  represent  propagation  through  a  turbulent  atmosphere.  The  index  of 
refraction  is  directly  coupled  to  variations  in  temperature  by  the  relation 
n  -  1  =-10"®  dT  [3].  In  this  approximation,  dT  is  an  incremental  change  in 
temperature  {°C)  within  the  isobaric  volume  element  of  air.  As  an  example,  a 
temperature  increase  within  the  isobaric  volume  element  of  dT~^°C  will 
decrease  the  air  density  in  the  region,  and  therefore,  decrease  the  index  of 
refraction  by  n-1=-10“®.  Conversely,  a  temperature  decrease  of 
dT  °C  will  increase  the  air  density  and  the  index  of  refraction  by 
n  -  1  =  1 0"®. 

The  index  of  refraction  n  is  related  to  a  dimensionless  spatial  increment 
dz  through  the  following  derivation.  By  definition  n  is  the  quotient  of  light 
speed  in  vacuum  and  light  speed  in  a  medium.  Specifically,  this  is  given  by; 
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n 


£ 

y 


(30) 


Differentiating  n  in  terms  of  v  and  regrouping  the  terms  gives 


dn  =  —^dv 


(31) 


By  approximating  that  v  is  nearly  equivalent  to  c.  Equation  (31)  simplifies  to 

dn  =  .  (32) 


This  difference  in  index  of  refraction  gives  the  light  speed  difference  dv,  and  the 
speed  difference  also  results  in  an  optical  phase  difference 


d(^ 


2ndvdz 

cX 


(33) 


over  the  random  element  of  path  dz.  Now,  recalling  that 
c/n  =  n  -  1  =  -10“®  dT,  then  the  index  of  refraction  can  be  incorporated  into  the 
propagation  code  by  the  approximation: 


d()) 


2k  X  10"® 


X 


dT 


(34) 


In  general,  beam  irradiance  deposited  on  a  target  is  dependent  on  the 
turbulence  encountered  along  the  propagation  path. 

A  phase  screen  operator  of  the  form  exp(/d(j)),  similar  to  the  lens  operator 
of  Section  C,  represents  the  effect  of  different  columns  of  volume  elements  with 
varied  indices  of  refraction.  Randomly  applying  phase  screen  operators  to 
segments  of  the  propagation  path,  results  in  the  simulation  of  a  random, 
turbulent  propagation  path.  Each  component  of  the  phase  screen  has  a 
random  t  vector.  The  amplitude  of  the  random  Jc  is  given  by  the  probability 
distribution  P{k)  «  /f"®'®[3].  An  exponential  distribution  of  characteristic 
wavenumbers  may  also  be  used  to  determine  the  Fourier  component  amplitude. 
The  direction  angle  %  is  used  to  determine  a  random  wavenumber  in  the  x  and 
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y  directions:  =  /fcos(27cx)  and  ky  -  ks\r\{2nx).  Each  phase  screen  is 

constructed  from  many  Fourier  components.  A  single  component  of  the 
numerical  phase  screen  operator  is  A(t)sin(/f;(.x  +  kyy  +  271x2)-  where  X2 
another  random  number  between  0  and  1.  Each  component  is  a  lens  with  a 
random  sinusoidal  surface  randomly  oriented  in  each  phase  screen  to  represent 
the  random  turbulence  over  the  propagation  path.  The  end  views  of  Figure  5 
and  6  illustrate  how  significantly  the  wavelength  affects  beam  scintillation.  A 
characteristic  wave  number  of  k^.  =  3  results  in  fewer  beam  spots  and  thus  a 
tighter  beam  concentration  on  the  target  then  Figure  6  which  has  kc=7  and  all 
other  parameters  equal.  This  is  in  general  agreement  with  the  experimental 
and  simulation  results  given  by  Gebhart  [9]. 


c=o. 02000 

nts512,  mc=64 
dE/E=:0.0000 


\ - i-i m 


|a|  (max) 


Figure  5.  Wavelength  affects  on  beam  scintillation. 


An  outcome  of  the  perturbations  along  the  propagation  path  is  that  the 
beam  phase  is  altered  and  its  energy  is  redirected.  Each  beam  spot  can  be 
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Figure  6.  Variation  of  characteristic  wave  number. 


approximately  the  size  of  the  transmitter’s  diffraction  limit  [7].  In  the  presence 
of  strong  turbulence,  the  beam  wave  breaks  up  into  groups  of  beam  spots. 
Thus,  the  phase  distortions  could  become  so  large  that  beam  focusing  would  be 
impossible.  This  is  of  concern  to  the  SFEL  model  since  it  is  crucial  to  have  the 
beam  focused  on  its  target.  A  comparison  of  Figures  2  and  6  reveals  how 
detrimental  turbulence  can  be  on  beam  quality.  In  Figure  2  the  beam 
propagated  through  a  homogeneous  atmosphere  and  focused  on  the  platform 
target.  The  simulation  in  Figure  6  maintained  the  same  parameters  but  added 
an  inhomogeneous  propagation  path.  The  nonhomogeneity,  or  turbulence,  is 
represented  by  the  total  phase  variable  =  600.  The  larger  the  value  of  the 
phase,  the  more  turbulence  encountered  by  the  propagating  beam. 
Additionally,  a  characteristic  wave  number  =7  was  added.  The  end  view 
reveals  a  dispersion  of  spots  concentrated  over  a  larger  area  when  compared 


to  the  end  view  of  Figure  2.  Thus,  the  turbulent  medium  makes  it  difficult  to 
focus  the  beam’s  energy  on  the  target. 

E.  THERMAL  BLOOMING  OVERVIEW 

The  effects  of  a  self-induced  phase  distortion  and  the  resulting  distortion 
of  the  laser  beam  irradiance  that  occur  when  a  laser  beam  propagates  through 
an  absorbing  medium  like  the  atmosphere  is  termed  "thermal  blooming"  [10]. 
The  absorbed  laser  beam  energy,  normally  a  small  fraction  of  the  total  laser 
beam  energy,  heats  the  medium,  causing  localized  gradients  in  the  density  of 
the  medium. 

Consider  a  SFEL  beam  propagating  through  the  atmosphere,  having 
some  absorption  at  the  laser  wavelength.  If  the  atmosphere  is  in  thermal 
equilibrium  and  the  SFEL  is  turned  on,  the  absorption  will  give  rise  to  local 
heating  which  produces  a  small  pressure  increase.  The  affected  atmosphere 
then  expands  at  the  speed  of  sound  so  as  to  restore  the  pressure  balance.  As 
the  medium  expands,  the  localized  density  will  decrease  with  a  corresponding 
drop  in  the  local  refractive  index.  If  the  cumulative  effect  of  the  negative 
refractive  index  change  is  large  enough,  it  will  result  in  a  change  in  the  shape  of 
the  laser  beam  profile.  Since  heating  usually  results  in  a  decrease  in  refractive 
index  of  the  medium  in  the  center  of  the  beam  where  the  heating  is  greatest, 
the  beam  is  usually  defocused  and  spread.  Thus,  the  beam  grows  in  size  or 
"blooms." 

Gebhart  [10]  enumerates  several  factors  that  influence  the  impact  thermal 
blooming  has  on  a  high  energy  laser  beam.  In  general,  these  include: 

1.  the  laser  beam  characteristics  (i.e.,  wavelength,  phase  and  irradiance 
distribution,  and  mode  such  as  continuous  wave,  single  pulse,  or 
repetitively  pulsed); 
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2.  the  kinetics  of  the  absorption  process,  which  determines  the  amount  of 
time  required  for  the  absorbed  energy  to  heat  the  atmosphere; 

3.  the  mode  of  heat  transfer  that  balances  the  absorbed  energy  (  i.e., 
thermal  conduction,  convection); 

4.  the  propagation  medium  and  scenario  characteristics  (i.e.,  path  length, 
optical  properties,  platform  speed,  slewing,  etc.). 

As  the  beam  passes  through  the  atmosphere,  there  may  be  locations 
along  its  path  which  may  allow  the  air  around  the  beam  to  be  stagnant  relative 
to  the  beam.  This  "stagnation  zone"  will  absorb  energy  from  the  beam  and 
heat  up.  The  heating  of  the  air  changes  the  index  of  refraction  of  the  air.  The 
energy  crossing  a  unit  area  will  be  proportional  to  the  square  of  the  field.  Since 
the  field  is  usually  stronger  in  the  center  of  the  beam,  the  air  in  the  center  of  the 
stagnation  zone  will  be  heated  more  than  the  edges.  The  effect  is  to  cause  the 
portion  of  the  beam  passing  through  the  center  of  the  zone  to  propagate  faster 
than  at  the  edges.  The  zone  acts  as  a  diverging  lens  defocusing  the  beam. 

As  an  example,  consider  a  1  liter  (10  cm)^  volume  element  of  air  with  an 
extinction  coefficient  of  0.05/km .  The  fraction  of  absorbed  power  in  the  volume 
element  is  0.05/km  x  10cm  =  5  x  10~®.  Thus,  if  the  laser  provided  1 
Megawatt,  the  absorbed  power  would  be  5  Watts  deposited  in  the  1  liter  volume 
element  of  air.  The  specific  heat  of  air  is  =  AE/{VAT)  ~  1J//°C,  where  AE 
is  the  amount  of  absorbed  energy,  V  is  the  volume  and  AT  is  the  temperature 
rise.  During  a  1  second  dwell  time,  the  temperature  rise  in  a  1  liter  volume 
element  would  be  5°C.  The  length  of  a  typical  stagnation  zone  is  50m  with  the 
center  of  the  laser  beam  passing  through  air  that  is  5°C  hotter  than  the  outside 
after  a  1  second  dwell  time.  The  path  length  difference  between  the  part  of  the 
beam  that  is  hot  and  that  which  has  not  is  calculated  in  the  same  way  as  (34). 
The  region  outside  of  the  laser  beam  is  not  heated  up  and  is  traveling  at  the 
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propagation  spoed  y  in  normal  air,  whoreas,  the  inside  of  the  laser  beam  is  hot, 
has  less  air,  and  is  traveling  at  speeds  greater  than  v.  Thus,  the  light  inside 
the  beam  travels  ahead  of  the  light  in  the  outer  regions  of  the  beam.  Thermal 
blooming  causes  a  phase  difference  between  the  inside  and  outside  of  the 
beam  in  proportion  to  the  laser  power  lal^.  The  phase  change  is  given  by  the 
following  relation  [3]: 

d^s  =  A(|)sla(x,y)|2  .  (35) 

The  simulation  program  calculates  the  effect  of  the  stagnation  zone  on  the 
real  and  imaginary  part  of  the  laser  field.  Figure  7  results  from  a  simulation  of  a 
stagnation  zone  at  t  =  0.6  with  an  effective  stagnation  zone  strength  of  Acfjg  =  2. 
The  stagnation  zone  strength  A<t>s  is  defined  similarly  to  the  lens  of  section  C 
however,  it  is  only  applied  to  the  stagnation  zone.  The  effect  of  the  zone  is  to 
cause  the  beam  to  defocus  and  is  not  random  in  Figure  7  because  no 
turbulence  is  included.  The  field  strength  in  the  center  of  the  beam  decreases 
relative  to  the  outer  parts  of  the  beam  cross-section. 

Figure  8  is  the  result  of  a  beam  propagating  through  an  atmosphere 
consisting  of  ns  =  100  turbulent  screens  with  a  characteristic  wavenumber, 
kc  =  4,  and  strength  A^^  =  600.  The  effect  of  the  turbulence  is  to  distort  the 
shape  of  the  beam  and  spread  it  over  a  larger  area  at  the  destination.  This 
Figure  combines  the  effect  of  turbulence  and  thermal  blooming  with  A(^s  =  2- 
The  distortion  of  the  beam  due  to  the  turbulence  is  apparent  as  well  as  the 
defocusing  effect  of  the  stagnation  zone.  The  result  is  less  concentration  of 
power  at  the  desired  point.  Since  blooming  is  affected  by  the  field  strength  of 
the  beam,  increasing  the  power  of  the  beam  will  only  exacerbate  the  blooming 
effect.  Therefore,  the  amount  of  power  focused  at  the  target  will  be  less  than 
expected  and  require  a  longer  time  for  the  beam  to  inflict  the  desired  affect  on 
the  target. 


23 


|side(0.99) 


top(0.99) 


end(0.52) 


C=0. 01000 
i  nt=1024,  nx=64 
dE/E=0. 0000 


0  X  1 0  .  I  a  I  (max) 

Figure  7.  Propagation  path  through  a  stagnation  zone. 
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Figure  8.  Propagation  path  through  a  stagnation  zone  and  turbulence. 
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Although  the  non-linear  blooming  process  can  dominate  the  beam  profile 
at  high  average  power,  linear  interactions  such  as  scattering  and  absorption 
combine  to  give  an  attenuation  or  extinction  coefficient  a  which  decreases 
beam  power  [3][10].  Scattering  due  to  aerosols  such  as  smoke,  mist,  and 
suspended  fine  particles  do  not  contribute  to  heating  of  the  air  in  thermal 
blooming.  The  primary  source  of  extinction  is  absorption.  The  absorption  of 
power  decays  as  P(z)  =  Poexp(-aT).  Where  Pq  is  the  original  power  level. 
Absorption,  clearly  the  basis  for  thermal  blooming,  is  very  sensitive  to 
wavelength.  A  SFEL  allows  wavelength  tuning  to  avoid  absorption  while 
maintaining  a  narrow  bandwidth.  This  provides  a  method  of  achieving  peak 
beam  irradiance  at  the  target.  The  wavelength  tunability  of  a  SFEL  is  an 
important  reason  for  selecting  it  over  a  conventional  laser. 
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III.  BOEING  1  KILOWATT  PEL 


Boeing  Corporation’s  Defense  and  Space  Group  has  designed  a  visible- 
to-infrared  wavelength  PEL  [13]  capable  of  delivering  1  kw  average  optical 
power  at  X,  =  0.6  pm  wavelength.  This  high-average  power  PEL  is  driven  by  a 
120  MeV  linear  accelerator  (linac)  generating  0.1  A  average  current  in  the 
macropulse  with  a  0.6%  duty  factor.  The  linac  uses  an  18  MeV,  433  MHz 
photo-injector,  and  a  1300  MHz  longitudinal  phase-space  linearizer.  The 
operational  characteristics  of  the  laser  are  presented  in  Table  1. 

The  first  operational  test  of  the  laser  is  expected  to  be  conducted  in  the 
visible  wavelength  region  during  fiscal  year  1996.  The  laser  and  supporting 
laboratories  are  housed  within  a  40,000  square  foot  PEL  research  complex  in 
Seattle,  Washington  [13].  This  facility  includes  a  12  MW  electric  power 
substation,  DC  power  vault,  accelerator  building  with  an  upper  radiation- 
shielded  floor,  and  general  purpose  laboratory  space.  The  electron  beam 
parameter  expectations  described  in  Table  1  result  from  many  experiments 
conducted  over  the  past  decade  in  support  of  the  Department  of  Defense 
Strategic  Defense  Initiative  and  the  Prench  PEL  program.  The  two  cavity 
photo-injector  has  been  operated  at  25  percent  duty  factor  and  32  mA  average 
current  [13].  This  demonstrated  capability  is  nearly  1000  times  that  performed 
by  any  other  photo-injector.  Additionally,  test  conducted  in  the  1980’s  at  the 
0.6  |xm  wavelength  indicate  that  the  expected  kilowatt  class  performance  is 
feasible.  An  extraction  efficiency  of  1%  with  strong  side-band  generation  was 
attained  using  a  much  more  inefficient  thermionic  injector  and  three  stage  pulse 
compression  [13]. 

Section  IV  highlights  the  theory  of  the  PEL.  In  Section  V,  applications  of 
the  Boeing  PEL  will  be  examined  with  an  emphasis  on  its  military  applications. 
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Parameter 

Value 

Energy  KE 

120  MeV 

Peak  Current  / 

755  A 

Undulator  Periods  N 

220 

Emittance  £„  {4*rms) 

30%  mm-mrad 

Undulator  Length  L 

5  m 

Energy  Spread  Sy/y 

0.004 

Undulator  Parameter  K 

1.31 

Beam  Radius 

0.016  cm 

Mirror  Separation  S 

20  m 

Rayleigh  Length  Zq 

250  cm 

Undulator  Wavelength  Xq 

2.18  cm 

Wavelength  X 

0.63  |j.m 

Pulse  Length  4 

7.0  ps/c 

Quality  Factor  Q 

90 

Table  1.  Characteristics  of  the  Boeing  1  KW  PEL  [13]. 


28 


IV.  BASIC  FEL  THEORY 


A  FEL  consists  of  three  basic  components:  a  relativistic  electron  beam,  a 
static  periodic  magnetic  field,  and  a  co-propagating  optical  field  [6].  A  stream  of 
relativistic  electron  bunches  pass  through  an  undulator  which  has  a  periodic 
magnetic  field.  Figure  9  is  a  simple  schematic  representation  of  a  FEL  and  its 
basic  components.  As  electrons  sinusoidally  traverse  the  undulator,  they  are 
accelerated  resulting  in  radiation.  Some  of  the  spontaneous  radiation  may  be 
stored  in  a  laser  resonator  formed  by  placing  two  curved  mirrors  at  the  ends  of 
the  undulator.  As  the  radiation  builds  up,  the  optical  field  in  the  undulator 
couples  with  the  "wiggling"  electron  beam  and  results  in  stimulated  emission 
and  coherent  radiation.  The  FEL  can  be  described  by  classical  electromagnetic 
theory.  This  results  from  a  high  density  of  photons  within  the  optical  wave  [12]. 
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Figure  9.  Schematic  representation  of  the  FEL.  From  Ref.  [11]. 
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A.  DIMENSIONLESS  PARAMETERS 


Recurring  combinations  of  physical  parameters  in  a  PEL  can  be 
summarized  by  a  few  dimensionless  parameters  [6].  This  discussion  will 
enumerate  many  of  the  dimensionless  variables  and  physical  processes  used 
throughout  the  following  sections.  Additional  parameters  will  be  explained  as 
needed  throughout  section  IV.  For  a  detailed  description  of  all  PEL  parameters 
see  reference  [6]. 

The  undulator  encompasses  a  series  of  magnets  arranged  to  permit  a 
periodic  magnetic  field  along  the  longitudinal  axis  of  the  undulator.  The 
dimensionless  undulator  parameter  \s  K  =  eBX(^{2Kmc^)  where  e  is  the 
electron  charge  magnitude,  B  is  the  undulator  root-mean-square  (rms)  magnetic 
field  strength,  ?io 's  the  undulator  wavelength,  m  is  the  electron  mass  and  c  is 
the  speed  of  light.  The  length  of  the  undulator  is  L  =  NXq,  where  N  is  the 
number  of  magnetic  field  periods  along  the  length  of  the  undulator.  The 
wiggling  electrons  interact  with  the  optical  field  along  the  entire  length  of  the 
undulator  for  a  dimensionless  duration  of  x  =  ct/L  where  t  =  0^1. 

The  relativistic  electron  energy  is  given  by  the  Lorentz  factor  y=  E/{mc^). 
Here  E  is  the  beam  energy  and  mc^  is  the  electron  rest  energy.  The  electron 
phase  ^  that  characterizes  the  electron  microscopic  motion  as  they  move  down 
the  undulator  is 

;  =  (/f  +  /Co)z(f)  -  (Of  ,  (36) 

where  the  optical  wavelength  wave  number  is  /c  =  2k/X,  the  undulator 
wavelength  wave  number  is  Xq  =  27t/Xo,  the  optical  wavelength  is  X,  and  the 
optical  frequency  is  &  =  kc.  The  only  dynamic  variables  in  Equation  (36)  are 
the  electron’s  longitudinal  position  z(f)  and  time  t.  In  the  relativistic  limit, 
k  »  Xq,  and  the  change  in  the  electron  phase  is  =  kAz  where  Az  is  a 
spatial  increment  along  the  longitudinal  axis  in  a  section  of  the  beam  one 
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optical  wavelength  long.  The  electron  phase  velocity  v(t)  is  the  dimensioniess 
time  derivative  of  the  eiectron  phase  as  given  by 

v{x)  ^  lix)  =  L[{k  +  ko)^,{t)  -  k]  •  (37) 


Where  the  electron  velocity  in  the  z  direction  is  given  by  cp^.  The  optical  and 
undulator  field  forces  are  resonant  when  v  =  0.  This  yields  the  resonant  optical 
wavelength  condition  which  demonstrates  the  wavelength  tunability  of  the  FEL. 
Specifically,  this  relationship  is  given  by: 


1  + 


(38) 


B.  ELECTRON  TRAJECTORY 

Electron  trajectories  in  the  undulator  (wiggler)  determine  the  character  of 
the  FEL  interaction.  It  is  the  electron’s  "wiggling"  motion  in  the  undulator  that 
provides  coupling  to  the  transverse  optical  field.  A  linear  undulator,  such  as  the 
"Thunder  Wiggler"  [13]  used  by  Boeing’s  1  kilowatt  average  power  FEL,  has  a 
series  of  magnets  along  its  length  which  provide  a  periodic  magnetic  field  and 
generates  linearly  polarized  light.  Another  type  of  wiggler  is  the  helical 
undulator  which  uses  coils  of  conducting  wire  wrapped  around  the  beam  line 
with  current  fiowing  in  opposite  directions  in  alternating  turns  [14].  This 
discussion  will  deal  mainly  with  the  linear  undulator  design, 

A  iinearly-polarized  magnetic  field  near  the  undulator  axis  can  be 

represented  by 

=  B[0,sin(/foZ),0]  ,  (39) 

where  B  is  the  magnitude  of  the  undulator  magnetic  field,  and  /Cq  =  2ii/Xo  is  the 
undulator  wave  number.  The  fields  within  a  linearly  polarized  undulator  in  the 
vertical  ^  direction  are  illustrated  in  Figure  10.  The  forces  acting  on  the 
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electrons  within  the  undulator  cavity  are  described  by  the  Lorentz  force 
equations: 


dt 


(1?) 


+  |x  +  fim) 


dy 

dt 


(40) 

(41) 


1 


(42) 


In  the  force  equations,  =  E[cos^,0,0]  is  optical  electric  field, 
=  E[0,sin'F,0]  is  the  optical  magnetic  field,  W  =  kz  -  (dt  +  ^  and  is  the 
electron’s  velocity.  Equation  (40)  represents  conservation  of  momentum  for  the 
system  and  is  commonly  referred  as  the  Lorentz  force  equation.  Equation  (41) 
describes  the  transfer  of  energy  and  (42)  is  the  Lorentz  factor.  Since  the  forces 
exerted  on  an  electron  by  the  magnetic  and  electric  components  of  the  optical 
field  nearly  cancel  each  other,  the  net  transverse  optical  force  is  significantly 
smaller  than  that  of  the  transverse  undulator  magnetic  field.  Thus,  the  magnetic 
field  in  the  undulator  is  the  dominant  field  moving  the  electrons  in  the 
transverse  direction.  Substituting  Equation  (39)  into  (40),  then  integrating  yields 
a  perfectly  sinusoidal  beam  trajectory  in  the  xz  plane  as  illustrated  in  Figure  10 
and  given  by  the  following  equation: 


[cos(;fo^).0.0]  •  (43) 

Y 

In  Equation  (43),  is  the  transverse  component  of  the  electron  velocity. 
Assuming  perfect  injection  into  the  undulator  cavity  makes  the  constant  of 
integration  zero.  Similar  approaches  can  be  used  to  describe  the  motion  of 
electrons  in  the  helical  undulator  [6]. 
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C.  RESONANCE  CONDITION 

Accelerated  electrons  in  the  undulator  of  the  PEL  radiate  energy.  These 
electrons  travel  along  the  longitudinal  axis  of  the  undulator  at  relativistic  speeds 
[6].  For  the  PEL  to  have  gain,  a  net  transfer  of  energy  from  the  electron  beam 
to  the  optical  field  must  occur.  This  energy  transfer  is  many  orders  of 
magnitude  greater  than  the  spontaneous  emission  contribution  per  pass  [16]. 
Therefore,  it  is  essential  that  stimulated  emission  occurs  preferentially  to 
stimulated  absorption.  For  an  optimum  exchange  of  energy  between  an 
electron  and  the  optical  field,  the  electron  must  oscillate  in  phase  with  the 
optical  field.  This  optimum  energy  exchange  occurs  when  the  electron  falls 
behind  the  optical  field  by  a  distance  of  one  optical  wavelength  within  the  span 
of  one  undulator  period. 

In  a  time  interval  Af,  the  optical  wave  will  travel  a  distance  cAf  and  the 
electron  will  travel  a  distance  p^cAf  where  the  axial  component  of  the 
electron’s  velocity  is  represented  by  p^.  Take  the  interval  At  to  be  the  time 
required  for  an  electron  to  travel  one  undulator  period  or  one  wavelength.  The 
light  travels  one  wave  length  X  ahead  of  the  slower  electrons  during  time  At,  so 
that 


X  =  {c-  p^c) 


(44) 


From  reference  [6]  we  know  that  p^  =  p|  +  Pf  and  recalling  the  Lorentz  factor 
from  Equation  (42),  the  axial  component  of  the  electron  velocity  can  be  found. 
Thus,  the  axial  component  of  the  electron’s  velocity  is 


1-pr- 


(45) 


Through  a  series  expansion.  Equation  (45)  can  be  approximated  as 
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(46) 


Pz 


+  7^Pi 

2f 


Squaring  the  magnitude  of  Equation  (43)  and  substituting  it  into  (46)  gives 

1  +K'2 


Pz  =  1- 


(47) 


Equation  (47)  relates  the  undulator  parameter  and  the  Lorentz  factor  to  the  axial 
component  of  the  electron’s  velocity.  From  this  simplification  we  can  reduce 
Equation  (44)  into  the  resonance  equation  which  is  Equation  (38).  The 
resonance  equation  reveals  the  capacity  of  the  FEL  to  be  tunable.  Since  the 
wavelength  of  the  optical  field  can  be  changed  by  adjusting  the  undulator 
parameter  and/or  the  energy  of  the  electron  beam,  the  FEL  can  be  continuously 
tuned  to  a  desired  wavelength  within  a  prescribed  bandwidth.  This  is  a 
significant  advantage  over  chemical  lasers  which  are  limited  to  a  design 
wavelength. 


D.  PENDULUM  EQUATION 

This  section  reviews  the  motion  of  the  electrons  resulting  from  the 
interaction  with  the  optical  field.  Electron  dynamics  in  the  FEL  can  be  modeled 
by  the  pendulum  equation  [6].  The  derivation  of  the  pendulum  equation  begins 
by  substituting  from  the  previous  section  into  (41).  This  yields 

T  =  ^^cos(^  +  (l))  .  (48) 

ymc 

Two  specific  cases  dealing  with  the  argument  of  the  cosine  function  in  (48)  are 
of  interest.  Firstly,  when  -k/2  <  C  +  <l>  <  7t/2  the  cosine  is  positive  and  the 
electrons  will  gain  energy  from  the  optical  field.  In  a  similar  fashion,  when 
ti/2  <  ^  +  ([)  <  3k/2,  the  cosine  is  negative  and  the  electrons  will  lose  energy 
to  the  optical  field.  Converting  the  time  increment  dt  in  (41)  into  a 
dimensionless  time  interval  dx  through  dt  =  Ldx/c,  Equation  (41)  becomes 


35 


(49) 


dy  _  eKLE 


cos(^  +  (})) 


dx  ymc‘ 
Differentiating  (47)  with  respect  to  x  gives 

(1  +  K^)J 


PzPz 


(50) 


where  the  notation  (..)  =  d{..)/dx.  Solving  for  and  substituting  in  (49)  gives 

h  =  +  .  (51) 


Differentiating  the  electron  phase  velocity  (37)  with  respect  to  x  gives 

OO  O 

^  =  L(/f  +  /fo)Pz  •  (52) 

Substituting  the  appropriate  relations  from  the  previous  sections  into  Equation 
(52),  and  assuming  /c  »  /Cq,  Pz  ~  "I  >  cind  imposing  the  resonance  condition  [6] 
results  in 

°°  ^  47cNe/CLEcos(^  +  ,]))  =  la  lcos(f;  +  (!>)  .  (53) 

'fmc^ 

In  the  above  equation,  the  dimensionless  optical  field  strength  is 
la  I  =  AnNeKLE/'fmc^.  Equation  (53)  is  the  electron  pendulum  equation  which 
governs  the  electron’s  phase-space  motion  within  the  undulator  and  under  the 
influence  of  the  optical  wave  [6]. 


E.  OPTICAL  WAVE  EQUATION 


By  adding  a  current  density  to  the  homogeneous  complex  optical  wave 
equation  described  in  Section  IIB,  the  complex  optical  wave  Equation  (1)  can  be 
rewritten  as  [6] 


2  dt\ 


W) 


(54) 


In  the  above  expression  the  transverse  electron  beam  current  is  .  Following 


36 


the  arguments  presented  in  Section  IIB  and  substituting  Equation  (4)  into  (54) 
yields 

=  -2)i7|  •  e,  .  (55) 

at 

and 

=  2TJ,-t2  ■  (56) 

ot 

Here  =  [cos'F -sin'P.O]  and  gg  =  [sinT.cos'F.O]  are  orthogonal  unit  vectors. 

The  transverse  electron  beam  current  is  the  sum  of  the  contributions  of 
each  individual  electron  [6].  Thus,  each  electron  contributes 

Tj,.  =  -  7) )  ,  (57) 

where  i)  is  the  position  of  the  ith  electron  and  6^®^(...)  is  the  three  dimensional 
Dirac-delta  function  [6].  Taking  the  dot  product  of  (57)  with  the  unit  vectors 
yields 


Assuming  a  constant  electron  density  p  and  summing  over  all  sample  electrons, 
the  electric  field  can  be  written  as 

E  =  -ZneKcp  <  >  >  (59) 

=  2neKcp  <  >  >  (®9) 

which  indicates  that  if  E  >  0,  the  laser  beam  energy  grows.  Specifically,  this 
will  happen  if  {^  +  ^)  ~  re.  Given  a  dimensionless  beam  current  [6]  of 
j  =  8N{enKLfp/i'fmc^),  and  the  dimensionless  optical  wave  amplitude  la  I 
from  Section  IVD,  Equations  (59)  and  (60)  will  yield 

la  I  =  -y  <  cos(C  +  (|))  >  ,  (61) 
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and 


I  =  -j-<  sin(^  +  (}))>  .  (62) 

la  I 

Equations  (61)  and  (62)  may  be  combined  to  yield 

a  =  -j  <  >  .  (63) 

Where  a  =  la  I  e'^.  Equation  (62)  reveals  that  an  increase  in  the  optical  field 
strength  will  result  in  a  decrease  in  the  optical  phase  evolution  rate,  while  an 
increase  in  the  dimensionless  electron  beam  current  will  result  in  an  increase  of 
the  optical  field  phase  evolution  rate.  Moreover,  in  the  absence  of  current, 
y  =  0,  the  FEL  optical  wave  envelope  does  not  evolve.  Similarly,  when  there  is 
current,  but  no  bunching,  and  thus  <..>  =  0,  there  is  also  no  change  in  the 
optical  field  [6]. 

F.  TRAPPED-PARTICLE  INSTABILITY 

At  the  onset  of  strong-field  saturation,  the  electrons  in  the  beam  tend  to 
travel  in  closed  phase-space  orbits  resulting  in  increased  FEL  efficiency. 
During  this  interaction  some  of  the  electrons  will  overtake  other  electrons  in  the 
beam.  This  results  in  diminished  bunching  and  gain.  A  direct  consequence  of 
this  is  that  some  electrons  will  become  trapped  in  deep  potential  wells  in 
phase-space.  This  phenomena  is  termed  the  "trapped-particle  instability."  It  is 
exhibited  when  either  the  FEL  amplifier  or  oscillator  evolve  to  power  levels 
beyond  the  onset  of  saturation.  Figure  1 1  depicts  a  group  of  electrons  traveling 
through  one  pass  of  the  undulator  and  becoming  trapped  in  the  potential  well  of 
phase-space. 

The  pendulum  equation  (53)  can  be  expanded  about  the  stable  fixed-point 
[6]  ^  =  n/2,  so  that  ^  =  n:/2  +  x,  where  the  spatial  increment  x<$iK.  A  solution 
describing  the  motion  of  the  trapped  electrons  is 
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Figure  1 1 .  Phase-space  evolution  showing  trapped-particle  instability  in  an 

untapered  undulator. 


CW  =  +  :^sin(v,t)  ,  (64) 

where  is  the  "synchrotron"  or  the  "trapped-particle"  frequency.  Under 

certain  conditions,  this  will  produce  sidebands  in  the  optical  power  with 
frequency  v  =  Vq  ±  and  wavelengths  shifted  away  from  the  fundamental  by 
AUX  =  Vs/{2nN)  [6]. 

The  interaction  of  electrons  along  the  length  of  the  undulator  can  be 
altered  by  modifying  the  properties  of  the  undulator.  Tapering  the  undulator  will 
alter  the  gain  spectrum  and  the  electron  beam  interaction  [6],  giving  the  FEL 
increased  power  at  saturation.  At  normal  saturation,  electrons  lose  sufficient 
energy  to  enable  them  to  drop  out  of  the  gain  spectrum  bandwidth.  As  the 
energy  of  the  beam  decreases  and  shifts  across  the  gain  bandwidth,  the 
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resulting  change  in  the  beam’s  phase  velocity  is  given  by  Av  =  AnN ~ -2%, 
where  Ay/y  is  the  change  in  the  beam’s  energy.  A  resonant  interaction  can  be 
re-established  by  tapering  the  undulator  field  in  order  to  extract  more  energy 
from  the  trapped  electrons.  When  operating  properly,  the  FEL  reaches  high- 
power  saturation  and  about  half  the  electrons  are  trapped  in  the  closed  orbits  of 
phase-space. 

The  dimensionless  synchrotron  frequency  for  the  trapped  electrons  in  a 
tapered  undulator  is  Vg  =  (lal^-  where  la  I  is  the  dimensionless  optical 

field  strength,  and  8  is  the  phase  acceleration  in  the  pendulum  equation  [6], 

oo 

^  =  la  lcos(C  +  <]))  •  (65) 

This  equation  shows  that  la  I  must  be  greater  than  5  in  order  to  trap  electrons 
[6].  In  an  untapered  undulator  Vg  =  VTal.  Comparing  Figures  11  and  12,  it  is 
apparent  in  Figure  12,  that  roughly  half  of  the  electrons  are  trapped  in  the 
potential  well  with  a  taper  6  =  47i.  Figure  12  illustrates  the  final  position  of  the 
electrons  in  the  separatrix. 

The  onset  of  the  trapped-particle  instability  can  be  seen  in  Figure  11.  For 
a  full  synchrotron  oscillation,  the  optical  field  strength  is  la  1=471^  =  40  [6], 
which  corresponds  to  a  separatrix  height  of  A\a\^'^  ~  25  [6].  This  causes  the 
electrons  to  travel  in  a  closed  orbit  about  the  fixed  point  (^,v)  =  (7i/2,0)  in  phase 
space  [6]. 

G.  DESYNCHRONISM  OVERVIEW 

This  section  only  presents  the  rudimentary  principles  of  desynchronism  d 
[6].  When  both  the  electron  and  the  optical  pulse  enter  the  undulator 
simultaneously  they  are  said  to  be  synchronized  when  d  =  0.  The  leading  edge 
of  the  optical  pulse  immediately  commences  overtaking  the  electron  pulse.  In 
weak  fields,  the  optical  wave  experiences  no  gain  at  the  beginning  of  the 


40 


Figure  12.  Phase-space  evolution  showing  trapped-particle  instability  in  a 

tapered  undulator. 

undulator  and  the  electron  bunching  is  slow  to  occur.  This  delay  in  the 
electron-optical  gain  mechanism  is  known  as  lethargy  [6].  The  initial  part  of  the 
optical  pulse  is  poorly  coupled  to  the  electron  beam  and  decays.  As  time 
progresses,  the  increased  gain  and  corresponding  bunching  of  electrons 
provides  the  necessary  interaction  for  field  coupling.  This  causes  the  trailing 
edge  of  the  optical  pulse  to  grow,  thereby,  distorting  the  optical  pulse.  The 
distortion  of  the  optical  pulse  continues  on  every  subsequent  pass  through  the 
undulator  as  gain  is  preferentially  deposited  on  the  trailing  edge  of  the  pulse  [6]. 
Essentially,  the  centroid  of  the  light  pulse  travels  slightly  slower  than  c.  Over 
many  passes  through  the  undulator,  gain  continues  to  grow  on  the  trailing  edge 
of  the  optical  pulse  resulting  in  eventual  power  decay.  This  occurs  when  the 
optical  pulse  centroid  moves  away  from  the  electron  pulse,  and  the  optical 


41 


wavefront  is  no  longer  driven  by  bunched  electrons.  Under  this  condition  the 
steady-state  power  goes  to  zero  [16]. 

Near  d  ~  0,  short  pulse  effects  restrict  coupling  and  the  steady-state 
power  is  zero.  Slightly  increasing  desynchronism  leads  to  a  sharp  rise  in 
optical  power,  where  the  electron  and  optical  pulse  repetition  rates  are  most 
nearly  matched.  Under  this  condition  the  FEL  could  become  subject  to  the 
trapped-particle  instability.  Increasing  d  further  causes  the  electron  and  optical 
pulses  to  stop  overlapping  over  many  passes  through  the  undulator.  Thus,  the 
interaction  of  the  electrons  with  the  optical  beam  deteriorates  resulting  in 
decreased  power  and  gain  [6].  In  summary,  too  much  or  too  little 
desynchronism  is  detrimental  to  FEL  power  and  gain. 
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V.  APPLICATION  OF  BOEING  FEL 


Applications  of  the  Boeing  1KW  FEL  can  be  categorized  into  three  groups: 
industrial,  research,  and  defense.  Industrial  applications  include  commercial 
isotope  production  and  pharmaceutical  material  processing.  Research 
opportunities  are  far  ranging  and  include  cancer  treatment  and  semiconductor 
development.  Defense  application  include  both  shipboard  defensive  and 
offensive  weapon  systems.  This  discussion  will  be  limited  to  the  use  of  the 
Boeing  FEL  as  an  Anti-Ship  Missile  Defense  (ASMD)  system. 

As  a  defensive  system,  the  FEL  is  an  excellent  ASMD  weapon  as  argued 
in  Section  I.  Assuming  Boeing  will  achieve  its  ambitious  objectives,  the  next 
concern  is  the  quality  of  the  optical  beam  exiting  the  FEL.  Boeing’s  high 
average  power  FEL  is  expected  to  yield  1kw  average  optical  power  at 
X  =  0.6  pm  wavelength.  This  section  will  focus  on  applying  the  principles 
presented  in  Sections  II,  IV  to  determine  if  the  necessary  beam  quality  can  be 
achieved.  The  focus  will  be  on  the  trapped-particle  instability  and  ways  of 
minimizing  its  effect  in  the  Boeing  FEL.  Boeing’s  plan  is  to  eliminate  the 
instability  in  order  to  narrow  the  power  spectrum  for  better  propagation  through 
a  transmissive  window  in  the  atmosphere. 

Figure  13  illustrates  the  trapped-particle  instability  in  an  untapered 
undulator  with  dimensionless  parameters  corresponding  to  the  values  of  Table 
1 .  The  dimensionless  current  density  j  =  658  is  derived  from  the  average 
current  [13].  The  parameters  ao  and  Vq  are  the  dimensionless  optical  field,  and 
electron  phase  velocity,  respectively.  Lastly,  Gq  is  the  induced  phase  velocity 
spread.  In  the  figure,  it  is  evident  that  most  of  the  particles  are  trapped  inside 
the  separatrix.  Thus,  the  necessary  interaction  between  the  optical  field  and 
the  electron  beam  is  not  taking  place.  In  fact,  the  electrons  are  absorbing 


43 


energy  from  the  optical  field.  The  peak-to-peak  separatrix  height  is 
~  4v5  =  4V5000  =  280.  Additionally,  the  synchrotron  oscillation  frequency 
is  Vs  =  70  giving  Vs/2%  ~  1 1  synchrotron  oscillations  along  the  undulator 
length  as  shown  in  Figure  13.  Oscillations  in  the  evolution  of  the  gain  G(t)  and 
the  optical  phase  (j)('T)  are  illustrated  at  the  right  side  of  the  figure  and  are 
caused  by  the  oscillations  of  the  electron  bunch  in  the  harmonic  closed  orbits 
around  the  stable  fixed  point  inside  the  separatrix. 


***  frt.  Phase  Space  Evolution  *** 

j=658  a  =5000  V^=10  N=220 

o  o 

Og=ll 

400 


V 


-400 

-7C/2  C  0  t  1 

Figure  13.  Phase-space  evolution  showing  trapped-particle  instability  in  an 

untapered  undulator. 

Figure  14  incorporates  a  10%  energy  taper  represented  by 
5  =  4%NAy/y~  88tc  [13].  The  energy  taper  dramatically  reduces  the  synchrotron 
oscillations  and  results  in  roughly  a  4%  increase  in  gain.  Additionally, 
approximately  half  of  the  electrons  are  trapped  in  the  potential  well.  These 
simulations  indicate  the  chosen  energy  taper  may  reduce  the  trapped-particle 
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instability.  However,  these  simulations  only  model  one  pass  through  the 
undulator.  Modeling  over  many  passes  through  the  undulator  is  necessary  to 
describe  the  performance  of  the  tapered  undulator. 


-7t/2  ^  37t/2  0  T  1 


Figure  14.  Phase-space  evolution  showing  trapped-particle  instability  in  a 

tapered  undulator. 

In  both  the  tapered  and  untapered  cases,  sidebands  can  form  around  the 
fundamental  frequency  at  Vo±nVs.  Figure  15  illustrates  the  results  of  a 
simulation  with  periodic  boundary  conditions  in  a  window  corresponding  to  one 
slippage  distance  along  the  beam.  The  dimensionless  current  is  j  ~  658  at  all 
points  in  the  window,  and  interacts  with  the  optical  field  for  n  =  400  passes. 
The  electrons  start  uniformly  spread  in  phase  along  z  with  a  phase  velocity 
spread  of  Cq  =  1 1 .  The  resonator  loss  on  each  pass  is  determined  by  Q  =  5 
corresponding  to  20%  loss  [13]. 
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Figure  15.  FEL  Wrap  Evolution  with  an  untapered  undulator. 

The  bottom-left  corner  of  Figure  15  shows  the  gain,  G{n),  at  the  end  of 
each  pass  through  the  undulator,  starting  at  resonance  vq  =  0  from  initial  optical 
field  ao  =  1,  with  noise  8a  =  1.  Other  modes  begin  to  grow  increasing  the  gain 
and  the  power  P{n)  shown  at  the  bottom  right.  As  the  electrons  slip  back  past 
an  optical  field  site  executing  synchrotron  oscillations,  they  continually  modulate 
the  wave  envelope.  The  middle-left  box  depicts  a  chaotic  optical  field  amplitude 
after  mode  saturation  and  the  growth  of  sidebands.  The  peak  field  amplitude  is 
la(z,n)l  =  2340.  The  middle-center  plot  shows  the  optical  power  spectrum 
P(v,n).  Power  P{n)  increases  dramatically  as  sidebands  appear  in  the  power 
spectrum  P(v,n).  The  final  power  spectrum  P(v)  shows  several  prominent 
sidebands  up  to  5v  =  200  away  from  the  fundamental.  The  frequency  spread  is 
large  at  ATuX  ~  bvl2nN  ~  ^5%  with  a  7%  shift  away  from  resonance.  The 
rectangular  tick  mark  points  at  the  center  of  the  final  power  spectrum. 
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The  electron  phase  velocity  spectrum  f{v,n)  at  the  top-right  corner  of 
Figure  15  shows  how  a  spread  is  caused  by  the  trapping  of  electrons  in  strong 
optical  fields.  The  rectangular  tick-mark  indicates  the  final  value  of  the  beam’s 
averaged  phase-velocities  and  shows  a  loss  of  energy  corresponding  to 
Ay/y  =  <Av>/47tA/  =  2%  with  a  large  spread  of  more  than  8%.  The  spread  in  the 
electron’s  phase-velocity  distribution  at  the  end  of  each  pass  f  {\,n),  increases 
as  the  strong  optical  fields  extract  more  energy  from  the  beam.  The  peak-to- 
peak  height  of  the  closed-orbit  region  of  phase  space,  4VTaT  =  4Vs,  gives  an 
estimate  of  the  induced  spread  at  each  pass.  In  normal,  single  mode  saturation 
the  electron  spread  is  only  Vg  =  2k,  after  the  trapped-particle  instability  occurs, 
the  peak  field  is,  la  I  =  2340  increasing  the  synchrotron  frequency  to  Vg  =  IStc. 

Figure  16  is  a  simulation  conducted  with  the  same  parameters  as  Figure 
15,  but  with  the  addition  of  a  10%  energy  taper  in  the  undulator  of  5  =  887c  and 
increased  Q  =  10  [13,6].  The  sidebands  are  not  eliminated.  Both  the 
sidebands  and  the  trapped-particle  instability  occur  causing  the  peak  field 
amplitude  to  be  chaotic  and  destroying  the  operation  of  the  tapered  undulator. 
The  average  wavelength  shift  is  Sk/X  ~  4%  with  a  wavelength  spread  of  =  8%. 
The  electron  beam  energy  spread  is  greater  than  10%  due  to  the  large  taper, 
but  most  electrons  remain  untrapped.  We  tried  simulations  starting  with  much 
stronger  optical  fields,  and  with  a  smaller  5%  taper  (8  =  44k),  but  could  not  find 
parameters  that  would  allow  a  single-mode  optical  field  so  that  the  tapered 
undulator  could  work  properly. 
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VI.  CONCLUSION 


The  development  of  a  high  average  power  FEL  for  military  applications, 
whether  shipboard  or  not,  represents  a  significant  advancement  in  technology 
over  present  weapon  systems  design.  The  FEL  has  significant  advantages 
over  conventional  kinetic  systems  and  other  classical  high-energy  laser 
systems.  The  rapid  response,  wavelength  tunability,  and  infinite  magazine 
make  the  FEL  a  highly  desirable  shipboard  weapon  system. 

The  proposed  Boeing  1kw  FEL  will  demonstrate  the  highest  average 
power  attained  to  date  in  an  FEL.  The  supporting  experiments  over  the  past 
decade  validate  their  high  expectations  of  achieving  high-average  power.  The 
simulations  of  Section  II  indicate  the  need  to  have  a  narrow  bandwidth 
spectrum  to  enable  propagation  through  a  turbulent  atmosphere  and  overcome 
the  affects  of  thermal  blooming.  By  tuning  the  wavelength  of  the  FEL  these 
atmospheric  losses  can  be  overcome  and  the  desired  energy  can  be  deposited 
on  target.  The  simulations  of  Section  V  indicate  that  the  Boeing  design  has 
potential;  however,  Boeing  must  correct  its  wide  bandwidth  problem. 
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